Tuning the Kondo effect with a mechanically controllable break junction 



J. J. Parks, 1 A. R. Champagne, 1 G. R. Hutchison, 2 S. Flores-Torres, 2 H. D. Abruha, 2 and D. C. Ralph 1 

1 Laboratory of Atomic and Solid State Physics, Cornell University, Ithaca, New York 14853 
2 Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853 

(Dated: February 6, 2008) 

We study electron transport through Ceo molecules in the Kondo regime using a mechanically 
controllable break junction. By varying the electrode spacing, we are able to change both the width 
and height of the Kondo resonance, indicating modification of the Kondo temperature and the 
relative strength of coupling to the two electrodes. The linear conductance as a function of T/Tk 
agrees with the scaling function expected for the spin- 1/2 Kondo problem. We are also able to tune 
finite-bias Kondo features which appear at the energy of the first Ceo intracage vibrational mode. 

PACS numbers: 72.15.Qm, 72.80.Rj, 73.63.-b 



The Kondo effect is a many-body phenomenon that 
can arise from the coupling between a localized spin and 
a sea of conduction electrons. At low temperature, a 
spin-singlet state may be formed from a localized spin- 
1/2 electron and the delocalized Fermi sea, leading to a 
correlated state reflected in transport as a zero-bias con- 
ductance anomaly pQ. This feature has been observed in 
devices containing lithographically-defined quantum dots 
[2] , carbon nanotubes [3] , and single molecules [4j [5j [6j [7] . 
Other more exotic Kondo effects involving higher spin 
states [H |9] , non-equilibrium effects [lOl [TTJ [12] , and or- 
bital degeneracies [13] have also been observed. 

While the Kondo effect has been intensely studied, 
previous experiments on nanotubes and molecules have 
lacked control over the molecule-electrode coupling. In 
this Letter, we address the influence of coupling on trans- 
port in the Kondo regime by tuning the Kondo effect 
in Ceo molecules with a mechanically controllable break 
junction (MCBJ) [14]. By varying the electrode spacing, 
we show that both the Kondo temperature Tk and the 
magnitude of the zero-bias conductance signal associated 
with the Kondo resonance are modified. These changes 
allow a determination of how the motion modifies the 
relative coupling of the molecule to the two electrodes. 
The normalized linear conductance exhibits scaling be- 
havior as expected within the theory of the Anderson 
Model. The same devices can also exhibit finite-bias 
inelastic Kondo features at an energy that corresponds 
to the lowest-energy intracage vibrational mode of Ceo- 
Changes in electrode spacing can tune the energy and 
amplitude of these signals, in ways that present a chal- 
lenge to existing theory. 

Our devices consist initially of continuous gold lines (32 
nm thick, 500 nm long, and with a minimum width of 50 
nm) suspended 40 nm above a 200-/im-thick Si wafer [15] 
(Fig. 1(a)). To incorporate Ceo molecules, we clean un- 
broken wires in an oxygen plasma (1 min at 30 mtorr, 
0.25 W/cm 2 ) to remove any organic contaminants. We 
then deposit 25 fiL of a 100 fiM solution of Ceo in toluene, 
wait 1 minute and blow dry. We cool the samples to 
1.6 K. The range of motion permitted by bending the Si 



substrate is generally too small to break the Au wires 
by mechanical motion alone, so we use electromigration 
to create a molecular-scale break in the wires before be- 
ginning studies as a function of mechanical motion [16]. 
We minimize the series resistance in our circuit during 
the electromigration process (the total series resistance 
is ~6 Q) and ramp a voltage stepwise until the resistance 
of a Au wire begins to change with time at constant V. 
This occurs consistently for V « 650 mV. We then hold 
V fixed until the wire breaks. After electromigration, we 
find that one or a few molecules can sometimes (~20% of 
194 devices studied) be found bridging the gap between 
electrodes, as inferred from the existence of a Coulomb 
blockade characteristic in the I-V curve rather than a lin- 
ear characteristic typical of bare junctions [4j[6]. In ap- 
proximately 3/4 of devices exhibiting Coulomb blockade, 
dl/dV displayed a peak at V = 0, which is a signature 
of the Kondo effect. We performed control experiments 
with more than 200 bare Au junctions as well as 76 junc- 
tions treated with toluene but without Ceo, and found 
a zero-bias feature in <2% of these devices: 3 bare Au 
junctions and 1 device exposed to toluene. 

We choose to study Ceo molecules because they are 
sufficiently durable to survive high temperatures present 
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FIG. 1: (color), (a) Scanning electron micrograph of a Au bridge 
suspended 40 nm above a Si substrate, (b) Schematic of the MCBJ 
apparatus, (c) dl/dV traces of Device A at various temperatures. 
Inset: A fit of the linear conductance to the NRG interpolation 
function (Eq. JlJ) yields a Kondo temperature of 28.2 ± 0.3 K. 
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during electromigration [17] and because previous work 
on single-molecule Ceo devices has observed the Kondo 
effect [6j 0. Photoemission studies of Ceo in contact 
with Au have shown that the molecule tends to gain an 
electron from the Au due to the molecule's high elec- 
tronegativity [THJ [19] , so that in our work Ceo is likely to 
often possess an unpaired spin in equilibrium, providing 
conditions needed for the Kondo effect. 

We can control the size of the inter-electrode gap in our 
devices by using a stepper motor attached to a pushing 
screw to bend the silicon substrate at cryogenic temper- 
atures (Fig. 1(b)). To calibrate the changes in electrode 
spacing, d, we measure the conductance of bare Au junc- 
tions as a function of motor turns and fit to the tunneling 
conductance expression G oc exp(— 2d^2m e (j)/h) 1 where 
<j) = 5.1 eV is the work function of Au. The mean calibra- 
tion over 14 bare junctions is 6.1 ±0.4 pm per motor turn; 
we apply this value to determine electrode displacements 
in identically-prepared devices containing Ceo- 

Figure 1(c) shows differential conductance curves of a 
Ceo device (Device A) at several temperatures. We ob- 
serve a prominent zero-bias peak that is suppressed with 
increasing temperature, in accordance with predictions 
for the Kondo effect. Shown in the inset of figure 1(c), 
we fit the linear response conductance as a function of 
temperature to an interpolation expression [20] that is a 
good fit to the numerical renormalization group (NRG) 
result [21] for conductance in the Kondo regime with the 
addition of a constant background conductance [22] , 
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using s = 0.22, and leaving Tk, Go, and G e \ as free pa- 
rameters. Our data are well described by this expres- 
sion, and we extract Tk = 28.2 ±0.3 K. The value for 
the Kondo temperature also agrees with that obtained 
by setting the full-width at half maximum (FWHM) of 
the base temperature zero-bias peak (5.14 ± 0.06 mV) 
to 2k B T K /e [3 , yielding T K « 30 K. 

Figure 2 displays the evolution of the zero-bias peak 
as we vary the electrode spacing for Device A (Fig. 2(a)) 
and for Device B (Fig. 2(b)). For each electrode spac- 
ing, we can determine the Kondo temperature from the 
FWHM of the zero-bias resonance, and we can also de- 
duce the relative coupling Tm/^l of the molecule to the 
more- and less-strongly-coupled electrode based on the 
magnitude of the linear conductance near T = 
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To describe the dependence on T, we adopt Eq. ([I]) and 
use f(T/T K ) = [1 ±T 2 /T K (2 1 / S - 1)]"" with s = 0.22. 
The value of the background conductance is determined 
from fits of the dl/dV as a function of V to a Lorentzian 
plus a constant. The low value of the peak conductance 
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FIG. 2: (color). dl/dV traces of (a) Device A and (b) Device B 
at various electrode spacings for T = 1.6 K. Insets to (a) and (b): 
dl/dV (in color scale) as a function of bias voltage and electrode 
spacing at T = 1.6 K. (c,d) The Kondo temperature and the rela- 
tive coupling Tm/IY as a function of electrode spacing for Devices 
A and B at T = 1.6 K. The uncertainty in the determination of the 
Kondo temperatures is smaller than ±0.4 K. 



of Device B indicates that the molecule is coupled quite 
asymmetrically to its electrodes. In figures 2(c) and 2(d), 
we plot the evolution of the Kondo temperature and rel- 
ative coupling as a function of electrode spacing for De- 
vices A and B. 

We find that as the electrode spacing is varied, Tm/IY 
increases by ~330% for a displacement of 0.55 A in the 
asymmetrically-coupled Device B, while the increase is 
just ^80% over a larger displacement of 1.1 A in the 
more-symmetrically coupled Device A. This suggests that 
as the electrodes are pulled apart, the molecule in Device 
B remains well-coupled to one of the leads so that the 
motion affects primarily IY, whereas in Device A, Tm 
and Tl are both modified, although not exactly equally. 
The background conductance G e i is always < 0.045 e 2 /h 
in Device A and < 6.4 x 10 -4 e 2 /h in Device B, and 
decreases by approximately a factor of 10 with 1.1 A 
motion in Device A, and by approximately a factor of 3 
over 0.55 A motion in Device B. 

In the more-symmetric Device A, the Kondo temper- 
ature can be tuned from ^30 K to ~23 K as the inter- 
electrode spacing is increased, whereas in asymmetric De- 
vice B, the Kondo temperature remains within a narrow 
range of 13-14 K. We can analyze these changes using 
the Haldane expression for the Kondo temperature in 
the limit of large charging energy U [23] , 
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where T = Tm + IY and So is the energy relative to 
the Fermi level of the localized state that produces the 
Kondo effect. We can expect both T and e$ to vary as 
a function of electrode spacing: V because the coupling 
of the molecule to at least one of the electrodes must de- 
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FIG. 3: (color), (a) Linear conductance of Device C as a function 
of temperature, fit to the NRG interpolation expression (Eq. JlJ). 
The extracted Kondo temperatures are (from the top trace to the 
bottom): 60.3 ± 2.4 K, 55.5 ± 0.9 K, 45.6 ± 1.9 K, 38.1 ± 1.2 K. 
Inset: dl/dV traces at xo + 0.7 A. (b) The normalized conductance, 
G(T/Tk)/Go, is a universal function of T/Tk- Scaled conductance 
data is compared to the NRG calculation (dotted line). 



crease as the electrodes are moved apart and Eq because 
break junctions generally exhibit large built-in electric 
fields even when V=0, so that motion of the electrodes 
produces a gating effect on energy levels in the molecule 
[T5] . If Tm and IY are significantly asymmetric, then 
varying the electrode spacing will likely have little effect 
on the overall T, as only the weaker coupling may change 
significantly. This regime applies to Device B, where the 
coupling ratio always exceeds 1500. The observation of 
only a small change in Tk as a function of electrode dis- 
placement for Device B is consistent with this picture. 
For Device A, we cannot distinguish the relative contri- 
butions of changes in T and So to the tuning of Tk, based 
on our data. In principle, a gate electrode that can inde- 
pendently adjust eo could help to disentangle the effects 
of adjusting electrode spacing. However, we find that 
the gate coupling for our device geometry is too weak to 
adjust eo measurably for devices in the Kondo regime. 

In figure 3(a), we show the temperature dependence of 
linear conductance for Device C at several different elec- 
trode spacings. In this device, the Kondo temperature 
could be tuned over a significant range, from 60.3 ± 2.4 
K (top curve) to 38.1 ± 1.2 K (bottom curve). In the 
Kondo regime, the conductance normalized to its zero- 
temperature value is expected to be dependent solely on 
T/Tk and thus to exhibit universal scaling behavior. In 
figure 3(b), we show that G(T/Tk)/Go does indeed ex- 
hibit a reasonable collapse onto a function that is in close 
agreement with the universal function predicted by NRG 
calculations [2T] . 

The evolution of transport properties as a function of 
electrode spacing is not always as smooth as we measured 
for Devices A, B, and C. In other devices, the conduc- 
tance could exhibit sudden changes, and zero-bias con- 
ductance resonances could fluctuate in and out of exis- 
tence. We ascribe these abrupt changes to motion of the 
Ceo molecule within the junction region. 

In addition to a zero-bias peak in dl/dV, in 5 out of 23 
devices with Kondo temperatures greater than 20 K we 



have also observed peaks in dl/dV at symmetric values 
of V near ±33 mV (specifically, at 29.6, 32.8, 33.5, 36.9, 
and 37.2 mV in the five devices); see figure 4(a). We 
did not observe any other similar features for \V\ < 60 
mV. The energy of 33 meV is known to correspond to 
the lowest intracage vibrational mode of isolated Ceo, 
in which the molecule oscillates between a sphere and a 
prolate ellipsoid shape (Fig. 4(a), inset) [24]. Previous 
investigations have shown that molecular vibrations can 
enhance dl/dV at energies corresponding to vibrational 
quanta [25j [26] [27]. For devices in the Kondo regime, 
coupling to a vibrational mode has been predicted to re- 
sult in an inelastic Kondo effect, which is manifested as 
sidebands in dl/dV at V ^ [28]. Finite-bias features in 
the Kondo regime have been observed in single-molecule 
transistors coupled to a vibrational mode [4j [6] [29], in 
quantum dots coupled to an applied microwave field [30 , 
and in quantum dots due to Kondo screening of excited 
states [flllll [311- 
Figures 4(b) and 4(c) show d 2 I/dV 2 for Devices D 
and E as a function of bias voltage and electrode spac- 
ing. Consistent with what we found for Devices A 
and B, as the electrodes are pulled apart the magni- 
tude of the zero-bias peak decreases more strongly in the 
less-symmetrically coupled (lower conductance) Device E 
than in the more-symmetrically coupled Device D. The 
strength of the sidebands is also different in the two de- 
vices; the satellite peaks are significantly more prominent 
in the more-symmetrically coupled Device D compared to 
the less-symmetrically coupled Device E, as predicted in 
Ref. [28 . As the electrode spacing, and hence the cou- 
pling asymmetry, is increased, the amplitude of the non- 
equilibrium peaks decreases in Device D, but the small 
peaks of Device E do not seem to be strongly modified. 
In both devices, the positions of the inelastic features 
increase in \V\ as the electrodes are pulled apart, sug- 
gesting that the mechanical motion increases the energy 
of the active vibrational mode. 

The changes in vibrational frequency as a function of 
mechanical motion are larger than what we anticipated 
based on molecular modeling. We performed calculations 
in Gaussian 03 using the PM3 semiempirical Hamilto- 
nian, which has been accurate in predicting the vibra- 
tional frequencies of fullerenes [32 . The geometries and 
harmonic frequencies were calculated under Cih symme- 
try using the native structure and by setting distance 
constraints on two atoms at opposite ends of the neu- 
trally charged Ceo cage to define the long axis from the 
native length of 7.092 A to 7.792 A. The calculations 
indicate that the five-fold degenerate H g (l) mode at 33 
meV is broken upon distortion into a set of two nearly- 
degenerate "short-axis" modes (involving motion perpen- 
dicular to the direction of stretching) which decrease in 
energy as the molecule is stretched and three "long-axis" 
modes which increase in energy (Fig. 4(d)). Only in- 
creases in energy are observed experimentally; we specu- 
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FIG. 4: (color), (a) dl/dV for Device D at T = 1.6 K showing 
satellite peaks near ±33 mV. Left inset: Schematic of the H g (l) 
intracage vibrational mode. Right inset: dl/dV for Device E at 
T = 1.6 K. (b, c) d 2 I/dV 2 as a function of bias voltage and elec- 
trode spacing for Devices D and E at T = 1.6 K. (d) Stretching de- 
pendence of mode energies from Devices D and E. Lines: short-axis 
(negative slope) and long-axis (positive slope) vibrational energies 
calculated using the PM3 semiempirical method. 

late that the long-axis modes may couple more strongly 
to electron transport. If we assume that the increase in 
molecular diameter is equal to the increase in the elec- 
trode spacing, so as to determine an upper limit for the 
estimated frequency shift, then we find that the mea- 
sured shift for Device E is comparable to the calculated 
frequencies of the long- axis modes, with the main devia- 
tions coming at large molecular diameter. However, the 
peak positions in the more-symmetrically coupled Device 
D shift more strongly than predicted, by roughly a fac- 
tor of 10. This discrepancy suggests that more rigorous 
theoretical work may be needed to understand electron- 
vibration coupling in single- molecule systems, when in- 
cluding coupling to the electrodes and the Kondo effect. 

In summary, we have demonstrated how the Kondo 
effect is modified by tuning the spacing between elec- 
trodes in mechanically controllable break junction de- 
vices containing Ceo molecules. We measure changes in 
both the Kondo temperature and zero-bias conductance 
that are in good accord with theoretical expectations and 
that allow us to characterize how the mechanical mo- 
tion changes the relative coupling of the molecule to the 
electrodes. We have also observed and tuned finite-bias 
Kondo features which appear at energies corresponding 
to an intracage vibrational mode of Ceo- We find that the 
vibrational energy can change more strongly as a func- 
tion of stretching than predicted by a simple semiempiri- 
cal Hamiltonian, thereby presenting a challenge for more 
accurate theory. 

We thank D. Goldhaber-Gordon, P. Brouwer, and A. 
Pasupathy for discussions, and K. Bolotin, F. Kuem- 



meth, and J. Grose for discussions and experimental help. 
JJP thanks the NSF for graduate fellowship support. 
This work was supported by the NSF (DMR-0244713, 
CHE-0403806, DMR-0605742), the NSF/NSEC program 
through the Cornell Center for Nanoscale Systems, and 
through use of the Cornell Nanoscale Facility /NNIN. 



[1] L. I. Glazman and M. E. Raikh, JETP Lett. 47, 452 
(1988); T. K. Ng and P. A. Lee, Phys. Rev. Lett. 61, 
1768 (1988). 

[2] D. Goldhaber-Gordon et al, Nature 391, 156 (1998); S. 

M. Crone nwett, T. H. Oosterkamp, and L. P. Kouwen- 

hoven, Science 281, 540 (1998). 
[3] J. Nygard, D. H. Cobden, and P E. Lindelof, Nature 

408, 342 (2000). 
[4] J. Park et al, Nature 417, 722 (2002). 
[5] W. Liang et al, Nature 417, 725 (2002). 
[6] L. H. Yu and D. Natelson, Nano. Lett. 4, 79 (2004). 
[7] A. N. Pasupathy et al, Science 306, 86 (2004). 
[8] S. Sasaki et al, Nature 405, 764 (2000). 
[9] W. G. van der Wiel et al, Phys. Rev. Lett. 88, 126803 

(2002). 

[10] S. De Franceschi et al, Phys. Rev. Lett. 89, 156801 
(2002). 

[11] A. Kogan et al, Phys. Rev. B 67, 113309 (2003). 

[12] J. Paaske et al, Nat. Phys. 2, 460 (2006). 

[13] P Jarillo-Herrerro et al, Nature 434, 484 (2005). 

[14] J. Moreland and J. W. Ekin, J. Appl. Phys. 58, 3888 

(1985); N. Agrai't, A. L. Yeyati, and J. M. van Ruiten- 

beek, Phys. Rep. 377, 81 (2003). 
[15] A. R. Champagne, A. N. Pasupathy, and D. C. Ralph, 

Nano. Lett. 5, 305 (2005). 
[16] H. Park et al, Appl. Phys. Lett. 75, 301 (1999). 
[17] M. L. Trouwborst, S. J. van der Molen, and B. J. van 

Wees, J. Appl. Phys. 99, 114316 (2006). 
[18] B. W. Hoogenboom et al, Phys. Rev. B 57, 11939 (1998). 
[19] C. T. Tzeng et al, Phys. Rev. B 61, 2263 (2000). 
[20] D. Goldhaber-Gordon et al Phys. Rev. Lett. 81, 5225 

(1998). 

[21] T. A. Costi, A. C. Hewson, and V. Zlatic, J. Phys. Con- 
dens. Matter 6, 2519 (1994). 

[22] M. Pustilnik and L. I. Glazman, Phys. Rev. Lett. 87, 
216601 (2001). 

[23] F. D. M. Haldane. Phys. Rev. Lett. 40, 416 (1978). 

[24] R. Heid, L. Pintschovius, and J. M. Godard, Phys. Rev. 
B 56, 5925 (1997). 

[25] H. Park et al, Nature 407, 57 (2000). 

[26] A. N. Pasupathy et al, Nano. Lett. 5, 203 (2005). 

[27] R. H. M. Smit et al, Nature 419, 906 (2002). 

[28] J. Paaske and K. Flensberg, Phys. Rev. Lett. 94, 176801 
(2005). 

[29] L. H. Yu et al, Phys. Rev. Lett. 93, 266802 (2004). 
[30] A. Kogan, S. Amasha, and M. A. Kastner, Science 304, 
1293 (2004). 

[31] A. Vidan et al, Phys. Rev. Lett. 96, 156802 (2006). 
[32] J. J. P. Stewart and M. B. Coolidge, J. Comput. Chem. 
12, 1157 (1991). 



